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Abstract 



We show the global existence of smooth solutions of a nonlinear partial dif- 
ferential equation modeling the dynamics of spinodal decomposition in diffusive 
materials. 
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1 Introduction 

In this paper, we consider an equation modeling phase separation in spinodal de- 
composition dynamics, which takes place in solid and liquid solutions under specific 
thermodynamical conditions. The initial stages of phase separation is revealed in the 
traditional Cahn-Hilliard theory jSj. However, Aifantis and Serrin in [3] suggested a 
generalization of the Cahn-Hilliard theory by including additional terms of interfacial 
stress. The derivation of their model equations is based on the balance laws of mass 
and momentum 



where u is the concentration, J is the flux of the diffusing material, F is the diffusive 
force and T is the symmetric stress tensor which includes the interfacial terms of a 
typical liquid-vapor phase transition. Combining the components of the tensor T in 
the one-dimensional case leads to the expression 
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where 5 is an interfacial coefficient and e is a short range deformity coefficient. Gen- 
erally, these coefficients are functions of concentration, but in this work, we assume 
them to be constants. The equation of state p(u) is assumed to be nonconvex with a 
single loop. On physical basis, the diffusive force F can be taken to be proportional 
to the flux J and its time rate, to incorporate inertia effects, that is 

F = M- l J + mJ t (3) 

where M is a mobility coefficient and m is a constant measuring the effect of inertia, 
which will be dropped in the current model, to obtain the relation 

J = -M(p(u) - eu 2 x - 8u xx ) x (4) 

Furthermore, we assume that the time dependent tensor T contains viscous relaxation 
terms of the form vu t . Combining the relaxation terms and letting M — 1 in P)l from 
the one-dimensional mass balance equation Q we arrive at the equation 

u t = p(u) xx + vu xxt - e{u 2 x ) xx - 5u xxxx (5) 

This is a fourth order nonlinear differential equation, solutions of which may explain 
the later stages of the transient spinodal decomposition process. Many aspects of 
the Cahn-Hilliard equation have been studied by Bates and Fife [3], Temam ^T] and 
Witelski [12J. The stationary and mechanical solutions of (0) have been investigated 
by Aifantis and Serrin in [21111] • In the current paper, we discuss the solvability of the 
equation (jSJ) for various boundary conditions, which are needed to explain the long 
time behavior and the evolution of spinodal decomposition process. 



2 Energy Estimates 

We consider the nonlinear evolution equation of spinodal decomposition for the den- 
sity u = u(x, t), that is 

u t = p(u) xx + vu xxt - e{u 2 x ) xx - 5u xxxx (6) 

on a bounded domain Q = [0,1] with the initial condition 

u(0,x) = u (x) e H 2 (n) (7) 

The equation (jHJ) is supplemented with either periodic boundary conditions (see 
Temam [9]): 

g(0,t) = ||(M) for i = 0,1, 2, 3 (8) 
or Neumann boundary conditions, see Bates and Fife [3]: 

^ (x ' t) = ^ (M) = forx = ' z - ( 9 ) 

In addition, we will assume that the equation of state p(u) E and grows linearly 
for \u\ > N for some large positive number iV and p'(u) changes its sign inside an 
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interval of displacement (phase separation). 



Local Existence and Uniqueness 

The derivation of local existence and uniqueness can be found for general differential 
equations in Henry [OJ, also presented briefly in Zheng Along the outline proof 
of local existence in time (see Henry [0]), we partition the differential operators into 
auxiliary linearized part and the remaining terms as follows: 



«* = ["(I " vdl)- l dl]{5u xx ) + [(1 - vdlY^lMu) + e{ul)) (10) 

We apply Fourier transform to the equation and treat the linear part as heat operator 
and denote the nonlinear part by / to obtain 

u = e~ 5 ^7(0 (11) 

We apply the contraction mapping theorem for sufficiently small time and energy 
estimation for the nonlinear part. The uniqueness can be shown by applying the 
above partition to the difference of two solutions w = u — v and concluding that w 
must be identically zero. 

However the crucial step in proving the global existence is to have a priori uniform 
estimates of the solution for any time T < +oo followed by the continuation argument 
(see Zheng [TU]). 

Global Existence 

Throughout this paper, ||.|| will denote the L 2 (Q) norm and c > will denote a 
generic constant that might depend on the initial data, e, 5, is, and possibly T but 
independent of t. The arguments in integrals will be omitted if they are clear. We 
prove the following theorem: Main Theorem The equation (Q|), equipped with the 
initial conditions |7|) and the boundary conditions (0) or (0J), has a unique global 
solution u £ C([0,T];H(Q)). 

Proof. The proof of this theorem is based on establishing a priori uniform estimates 
on the solution u. We group these estimates into two lemmas. 



Lemma 2.1 For any t £ [0, T] we have 



rT 

sup (|M| 2 + |K|| 2 ) + / \\u xx fdt<c. (12) 

0<t<T 



Proof. Multiply equation (JBJ) by u and integrate by parts over Q to obtain: 

pi pi pi pi pi 

I uu t dx = / up xx dx — vl uu txx dx — e I u{u 2 x ) xx dx — 51 uu xxxx dx (13) 
Jo Jo Jo Jo Jo 

We evaluate each term using the boundary conditions (jHJ) or Q and the restrictions 
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on the pressure p, as follows 
fi 

UUtxxdx = uu tx 
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+ / u x u( x ) x dx 
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< k u x dx = k\\u a 
o 



We combine these estimates to obtain 

1 



u\\ + u\\u x \\ ) t + S\\u xx \\ < k\\u a 



(14) 



(15) 



Integrate the inequality (fTB^) over time interval [0,T] and invoke Gronwall's lemma 
to deduce the required estimate (fl^jl . 



Lemma 2.2 There holds 



sup (\\u t \ 

0<t<T 



\u tx 



\u txx \\dt < c. 



(16) 



Proof. Differentiate © with respect to t and multiply by u t to get 

u t u tt = u t p xxt + vu t u xxtt - eu t {u 2 x ) xxt - 5u t u xxxxt . (17) 

Integrating (fTTj) by parts over Q taking into account the boundary conditions and the 
estimates from lemma [2~T1 yields the following relations 

u t u xxtt dx = (u t u xU )\ l - / u xt u xtt dx = 

Jo 

1 

'2' 



u 2 xt dx = -~|Kt|| 2 . 



o 



-J u t (u x ) xxt dx = (u t (u x ) xt )\ + j u xt (u x ) xt dx 

= - u xxt {u x ) t dx = 
Jo 

— — 2 / u XX fU x u X fdx . 
Jo 
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Apply the Cauchy-Schwartz inequality to the last integral to get 

u t {ul) xxt dx < 2x I u 2 xxt + - / u 2 x u 2 xt dx 



o X Jo 



< ZxhxxtW + -\\u xt \\U / u 2 dx (19) 
X Jo 
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< 2x\\u xxt \\ 2 + —\\u xt \\ 2 L oo \ u 2 x dx. 

X Jo 

We evaluate the last term in (flHj) by applying the Young inequality 



\u x t\\h^ < ciiiM^tH 1 / 2 !!^!! 1 / 2 + c 2 ||w x t|| 

< CiX 2 ||M^|l + "AlKJ + Ca||"srt||- (20) 



Regrouping these estimates to obtain 

ll^tllloo < cx 2 \\u xxt \\ 2 + c(x)\\u xt \\ 2 . (21) 
where c(x) is a function of x- Similarly, we evaluate the remaining terms to obtain 

1 k 
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u t p(u) xxt dx 



< -\\u t \\ 2 + k X \\u xxt \\ 2 (22) 
X 



and 



UtU XXX xtdx — 1 1 1ia;a;t 1 1 (^3) 







Substituting these estimates and selecting a small enough x we arrive at the following 
inequality 

-(||u t || 2 + v\\u xt \\ 2 ) t + ^oIK^H 2 < Cill^H 2 + c 2 ||w x .J 2 , (24) 

where is a positive constant. Integrating ()24|) over [0, T] and invoking the Gron- 
wall's lemma. We conclude (jTHJ). 

Remark 1 We can derive additional energy estimates for higher order derivatives of 
the solution of equation (0) by interpolation relations and requiring suitable degree of 
smoothness of initial data. 
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